Analytical solution of a generalized Penna model 
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In 1995 T. J. Penna introduced a simple model of biological aging. A modified Penna model 
has been demonstrated to exhibit behaviour of real-life systems including catastrophic senescence 
in salmon and a mortality plateau at advanced ages. We present a general steady-state, analytic 
solution to the Penna model, able to deal with arbitrary birth and survivability functions. This so- 
lution is employed to solve standard variant Penna models studied by simulation. Different Verhulst 
factors regulating both the birth rate and external death rate are considered. 

PACS numbers: 87.23.-n, 87.10. +c 



I. INTRODUCTION 

While nothing in life is certain except death and taxes 
0, only death is universal. Death and the preceeding 
period of functional decline are a fate to be endured by 
all from medfLies to men 0, OJ 0, Q • The steady decline 
in functional ability of an organism over time is known as 
senescence or aging. The phenomenon of aging is of ob- 
vious interest and has attracted attention from biologists 
and physicists for some time Q- 

The establishment and maintaintenance of harmful be- 
haviour by natural selection would appear to defy expla- 
nation. Medawar proposed 0, @ that the strength of 
selection on survival related genes is dependent on the 
age at which the genes exert their effects. For genes that 
express themselves late in the life of an organism there is 
less impact on the dwindling population than for a gene 
whose effect is expressed earlier in life. 

It has been suggested % that favourable mutations 
that act to increase survivability can be used to account 
for senescence. Such mutations will, under natural selec- 
tion, increase the survivability at early ages converting 
an initially constant mortality rate into an age depen- 
dent one. It seems improbable that this mechanism can 
provide a full explanation of aging as positive mutations 
are so very rare compared to harmful ones. Instead a bet- 
ter understanding can be gained by considering processes 
through which genes with harmful effects are introduced. 

There are two such theories: antagonistic pleiotropy 
and mutation accumulation [Tol ITU - According to an- 
tagonistic pleitropy, senescence occurs as a result of mu- 
tations that increase the functional ability of the young 
and decrease that of the old. Mutation accumulation pro- 
poses that aging occurs due to mutations that are initially 
harmless but take effect at later stages in the life of an 
organism. In either case, the force of natural selection is 
reduced once an organism ages beyond its point of repro- 
ductive maturity so the effects of either will be confined 
to older ages. 

The Penna model is the most commonly used model 
of aging through mutation accumulation |l2l and is 
ideally suited to computational implementation. Using 
Monte Carlo methods |14|, the model predicts features 
found in real ecological systems such as the catastrophic 



senescence of Pacific salmon [l^. Analytical work by 
Almeida et al 16] on a theoretical approach to biolog- 
ical aging can be adapted to apply to specific cases of 
the Penna model. This work presents a solution to a 
generalized Penna model, in particular one that allows 
incorporation of arbitrary survivability and birth func- 
tions. Subtle modification to the survivability function 
has been found |l8j to demonstrate a mortality plateau 
at older ages, a result that had so far eluded theories of 
mutation accumulation. 

To ensure that the population is in a steady state, 
the total population is controlled through the use of a 
Verhulst factor |l9|, |20|. Traditionally this has been a 
genome-independent chance of death for every individ- 
ual regardless of age. This model also considers a birth 
rate that decreases as the population grows and resources 
become more scarce |2l| . In agreement with earlier work 
|16| . we find a maximum permitted genetic lifespan and 
predict the existence of catastrophic senescence for or- 
ganisms whose reproductive life is terminated by an up- 
per age limit. The model is extended to a continuum 
case, which is explicitly solved. 



II. THE PENNA MODEL 

The Penna model as formulated by Penna represents a 
genome by a single string of Is and 0s. Time is treated as 
a discrete variable. A 1 on a site i along the string means 
that the organism develops a disease at age i. Once an 
organism develops a number, T, of diseases it dies. At 
each time step an organism reproduces with probability 
b. The offspring's genome is a copy of its parent's with 
a probability m of each bit mutating into a 1. Positive 
mutations are rare in nature and ignored in this model, 
there is no possibility of a 1 mutating into a 0. The bit 
string is traditionally 32 bits long for ease of computa- 
tional implementation. The finite-length bit string is an 
artifact of simulation and is discounted in our analysis 
where there is no need for such a restriction. Along simi- 
lar lines to Almeida and Thomas |17| we consider a fixed 
probability of mutation occuring on any site. The bit- 
string sites are labelled so that there is a zeroth site that 
is read as soon as the organism is born. 
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A. The Solution of a simple Penna model 

Consider a simple Penna model where an organism dies 
after a single disease (T = 1) and can reproduce with 
equal probability at any point during its life. An indi- 
vidual organism can be characterized uniquely by two 
variables, its age x and its string length I. The string 
length is the location of the 1-bit on the string and cor- 
responds to the number of time steps for which the organ- 
ism lives. To produce an organism with string length I 
either a perfect copy of a length I organism or a mutated 
copy of an organism with longer string length must be 
born. To produce an organisms of string length 1, an or- 
ganism must give birth, with probability b, and the first 
I sites on the offspring must go unmutated, each with 
probability (1 — to). For mutated offspring, the parent 
string must be longer, the parent must give birth, with 
probability 6, the first I sites on the offspring must go 
unmutated and one site must be mutated with proba- 
bility m. As all organisms are capable of reproduction, 
mutated and perfect copies of organisms of any age must 
be taken into account. 

In our notation rij(x,l) is the number of organisms 
with age x and string length I at time-step j. We define 
e~P to be (1 — to). New organisms are produced as mu- 
tated or unmutated copies of organisms in the previous 
time step. 

oc 

rc i+ i(0,l) = be- pi J2 n A x > 1 ) (!) 

oc oc 

+ mbe- ffl J2J2 n i( x > 1 '^ 

V>1 x=0 

For a steady-state there is no difference between n(x, I) 
at different time steps so the time step indices can be 
dropped. The simple Penna model is constructed in such 
a way that an organism with string length I lives for I time 
steps. The probability of giving birth during each of these 
time steps is a constant b. Thus the sum over all ages is 
a sum from to I of n(x, I) which in the steady-state can 
be written as I x n(0, 1). Defining n(l) — I x n(0, 1) the 
equation now reads 



= be-P l n(l) - ^ 



ibe 



-pi ' 



v>i 



n{l'). 



(2) 



A similar equation can be written for n(l + 1) manipu- 
lation of both will eliminate the sum over I' and give a 
recursion relation 



n(l + 1) _ I + 1 



i0l 



bl 



I eW+i) - 6(1 + l)e- 



(3) 



If this expression is to be usefully employed the steady- 
state interdependence of b and (5 must be examined. In 
the steady-state it is required that, in the statistical limit, 
the population size and distribution remain unchanged 
over time for constant values of b and j3. 
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FIG. 1: Lifespan distribution for Z max = 30 Analytical results 
(x) are compared with those from simulation (□). Simulation 
size 10 7 . 



As the string length I gets longer, the probability of 
an unmutated copy being produced is reduced by an ex- 
ponential factor e~@ l . For an organism to be able to 
produce a single perfect copy of itself during its lifespan 
lbe~@ l — 1 If any sub-population is capable of maintain- 
ing itself, there must be no contribution to that sub- 
population from mutations. As all sub-populations will 
have contributions from mutation from longer strings, 
it follows that a sub-population capable of maintaining 
itself must have the longest string length in the popula- 
tion, Z max . For all other sub-populations with I < l max 
the probability of an organism producing a perfect copy 
of itself must be less than one to avoid population ex- 
plosion when contributions from mutation are taken into 
account. 

It is sufficient to state that 



(1, 



l)6e 



-0(h 



-i) 



< 1. 



(4) 
(5) 



The constraints on l max lead to a maximum sustainable 
value and a corresponding steady-state birth rate 



max ^ 
b = 



1 



1-e-f 
1 



(6) 
(7) 



In agreement with [l6| we have predicted the existence 
of a maximum sustainable genetic lifespan, see Figs. 1-6. 

The system has a range of possible values for Z max and 
will adopt one depending on the dynamics and the initial 
state of the system. In simulation, a population is often 
regulated through the use of a Verhulst factor |T3 . fT^. l20| . 
This ensures that if a population will find a steady-state 
configuration. Such simulations also regulate the total 
population which is realistic in any system with finite re- 
sources. It has been suggested [21| that by regulating the 
birth rate, the system will adopt a biologically realistic 
equlibrium. The Simple Penna model solved here can 
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be used to explain behaviour of these simulations that 
replace the constant birth rate with a population depen- 
dent one. Where N is the total population, A a constant 
of the simulation and N max is the maximum population 
the simulation will allow, the birth rate at any time-step 
i is given by 



h = A 1- 



Ni 



N„ 



(8) 



In the steady-state the total population will be that re- 
quired to set this dynamic Verhulst birth rate to the value 
required by the earlier equilibrium conditions. 



N = N n 



1 



AL 



(9) 



B. 



The solution of a Penna model with 
reproductive threshold ages 



Frequently the Penna model used [F| has birth cut- 
offs so that an organism begins giving birth at an age, Ib 
up to an age, Is- When able to reproduce, an organisms 
does so at rate b as before. The equation for step-wise 
evolution of an age zero sub-population remains the same 
except that the birth rate is now a function of the organ- 
isms age. An organism is no longer capable of reproduc- 
ing in every time step, but only those for which its age 
is greater than or equal to Ib and less than Is- The sum 
over ages must now sum over the age-dependent birth 
term as well: 



n(0,l) 



-pi 



(10) 



rn e 



x=0 

CO oo 

-pi 



^2^2b x n{x,l'). 



l'>l x=0 



As in the simple Penna model n(l) is defined to be 
I x n(0, 1) as the survivability function is unaltered. Sum- 
ming over reproductive ages leads to a simplified steady- 
state equation 



= be-P l x(l)n(l) 



n(i) 



mbe- 0l J2x(l')n(l'), (11) 



v>i 



where %(/) is defined as ^"^(o")^'^ an d is given by 



X(0= for l<l B 

X(l) = I-Ib for l B <l<ls 
x(l) = Is - Ib for I > l S - 



(12) 



Employing the same method as was used to solve the 
simple Penna model, a recursion relation can be obtained. 



n(l + 1) _ I + 1 



bx(!) 



I e P(i+i) _ b x (i + i) e -/3 



(13) 



0.08 r- 

0.07 - 

c 0.06 - 
o 

aj 

3 0.05 - 

Q. 

i 0.04 - 
o 

■1 0.03 - 

o 

nJ 

£ 0.02 - 
0.01 - 




10 15 20 

Genetic Lifespan L 



25 



30 



FIG. 2: Lifespan distribution for i m ax = 30, Ib = 5. Analyti- 
cal results (x) are compared with those from simulation (□). 
Simulation size 10 7 . 



In the case of no birth thresholds, = I and we retrieve 
©■ 

The steady-state conditions for a population will have 
been altered by the introduction of birth threshold ages. 
The conditions for Z max are, by the same reasoning as 
before: 



{In 



l)be 



''max 

0(1 



be-^x(ln 

max 



i) 



= l 
< l. 



(14) 
(15) 



There are two distinct non-trivial cases for the position of 
Z max , either it is below the upper reproductive boundary 
or is contained within the reproductive window. 

Consider Z max < Is- The boundary conditions can be 
written as: 



(Zmax " l B )be-^ 

After rearrangement, this yields 



= 1 
< 1. 



< 



1 + {e-P - 1)1 



with a corresponding birth rate of 

„/3/ ma x 



b = 



B 



(16) 
(17) 



(18) 



(19) 



Consider now the case of i max 
conditions become 



b(l s - l B )e' pl — = 1, 
b{l s ~l B )e- p{l — - 1] < 1. 



> ^5. The boundary 



(20) 
(21) 



Requiring e 13 < 1. This cannot be satisfied as (3 is a muta- 
tion rate and consequently a positive number. As a result 
^max is not permitted to be greater than the second birth 
threshold Is for any value of mutation rate or birth rate. 
A steady-state population will not maintain organisms 
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that are no longer able to contribute to reproduction. 
This result does not apply to organisms who have some 
social structure and assist in rearing young once they 
themselves are no longer of a reproductive age. Such co- 
operative complications are not taken into account 
here. 



C. The solution of a Penna model with an external 
death rate 

Consider the simple Penna model. An external death 
rate is introduced into this model so that in any time 
step an organisms has a chance 7 of dying independent 
of its bit-string composition. We define the survival rate 
a as 1 — 7. Counting contributions from unmutated and 
mutated reproduction, the steady state equation is 

00 00 00 

n(0, = be' 131 J2 n ( x > + mbe- 131 J2 '')• ^ 



x=0 



!'>/ x=0 



The sums over all ages can be evaluated by taking into 
account the external source of death. 



1-1 



n(0,0- 



(23) 



x=0 



x=0 

I -a 1 



7 



Defining n(l) = 
be re-arranged to give 

= be~ m n(l) 



-n(0,l) the steady-state equation can 



1 - a 1 



i{l) + mbe- 0l Y^n{l'). (24) 



v>i 



The same approach is used as before to give an recursive 
relationship between nil) and nil + 1). 



n(l + 1) 
n{l) 



1 — (7 



1 + 1 



I -a 1 

ye? 1 -b{l + a l ) 
7e /J(J+i) - 6(1 - tr'+^e-/ 3 



(25) 



In the limit of a vanishing external death rate, (small 
7) power-series expansion of these expressions will, to 
leading order, give the conditions from the simple Penna 
model. ©. 

To find the steady-state relationship between b, (3 and 
7 the same conditions are imposed as before. n(Z max ) 
must be self sustaining and all lower string lengths are 
partly reliant on mutation. 



= 1 



(26) 
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FIG. 3: Lifespan distribution for Z ma x = 30, 7 = 0.02. Ana- 
lytical results (x) are compared with those from simulation 
(□). Simulation size 10 s . 



These conditions give the birth rate and limit the value 



of L 



< 



lner 



ye h 



l-(T l 



(27) 
(28) 



In the limit of a small external death rate these expres- 
sions become those from the simple Penna model case 

An external death rate is commonly used in simula- 
tions 0, ^| as a Verhulst factor to regulate the popu- 
lation. In any time step an organism has probability of 
death given by J** The preceeding analysis can be used 
to explain how systems using this Verhulst factor behave. 
In the steady-state the external death rate provided by 
this Verhulst factor must satisfy the above conditions for 
stability at and around ( max For specified 6, (3 and /max, 7 
can be found numerically and the total population found 
from N = 7A r max - 



D. The Solution of a Penna model with an 
External Death Rate and Reproductive Threshold 
Ages 



It is a relatively trivial matter to extend the analysis 
of a Penna model with an external death rate to incor- 
porate birth threshold ages. This type of model has been 
extensively studied in simulation so is of sufficient inter- 
est to consider separately. Only the lower birth threshold 
needs to be considered as an upper threshold, as demon- 
strated, will only serve to artificially lower the maximum 
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FIG. 4: Lifespan distribution for Z max = 30, Ib = 5, 7 = 
0.02. Analytical results (x) are compared with those from 
simulation (□). Simulation size 10 8 



allowable string length, Z max - As ever, in the steady-state 

00 

n(0,l) = e -/ " £ Mfo (29) 

oc 00 

+ me^ l3l ^2^2b x n(x,r). 



V>1 x=0 



Taking external deaths into account when summing over 
reproductive ages we define x 7 (Z) as 



X 7 (0 = for I < l B 



(30) 



n(l) can be given 

n(/ + l) _ l-a l+1 e^- 



with n(£) given by - — ^^Ln(0, Z) a recursion relation for 



(31) 



n(l) I -a 1 e^ l+1 ) ~ b Xl {l + l)( 



E. The Solution of a Multiple Disease Penna model 

An organism in the simple Penna model dies upon en- 
countering a single 1 on its bit-string, models are fre- 
quently set up so that an organism must develop T dis- 
eases before death. Each site on the string is occupied 
by one or no deleterious mutations, multiple mutations 
on a single site are not allowed. 

The relevant part of an individuals string is that con- 
taining the first T Is, as bits after this point are irrele- 
vant. Rather than uniquely classifying an organism by its 
genetic lifespan I as was done in the single mutation case 
we must now specify the position of each of the first T 
mutations on the organism's bit string. The only other 
property an organism has is its age x. Any individual 
can thus be uniquely classified by its type x; ii, 1% . . . It 



where n(x; h, 1% . . . It) is the number of the specified or- 
ganisms. The position of the final disease on the string 
determines an organism's age of death, the position of 
the other bits play no part in this, nor do they deter- 
mine birth rate. This inspires the ansatz that n(x; {I}) 
has no dependency on the positions of the non-terminal 
diseases. 

Contributions to a child sub-population from a sin- 
gle mutation come from all organisms with T — 1 bits 
in common with the child sub-population, and the final 
bit at a position V > It- For T = 4, where l±, 1%, I3 
and I4 are the deleterious bit positions in the child sub- 
population, there will be contributions from mutation 
from n(x;h,l2,l 3 ,l'), n(x;lx,l 2 ,k,l'), n(x;h,l 3 ,l4,,l') 
and n(x;l2-,l3,l>i-,l')- Our ansatz means that all these 
terms are the same as they are dependent on only the 
age x and the position of the terminal bit I'. Where 
nr(x, I) is the number of organisms with age x and ter- 
minal mutation at site I the number of organisms capable 
of contributing through a single mutation is TriT{x,l'). 
In a similar vein to the single disease Penna model we can 
now account for all birth terms and generate a steady- 
state equation. 



n T (0,Z) = beW-T+V^nrfal) (32) 

x=0 

00 00 

+ Ttnie^- T+1 )j;£nr(s,!') 



V>1 x=0 



This can be solved in the same manner as that for the 
simple Penna model, where nr(l) — I X nr(0, 1) 



n T (l + l) 
n T (l) 



l + l 



(33) 



-M-T+l) _ bl 



g-W + l-T+l) _ b{ J + 1)(e /3 + T - 1) ' 

UtQ) is not of direct interest in determining age dis- 
tributions or mortality rates as it is only one of sev- 
eral configurations of non-terminal mutations. This is 
amended by summing over all possible configurations so 
that n(l) = C t _ x tit{1)- 

The steady-state correspondence between 6, (3 and Z max 
for arbitrary T is 



max ^ 
b = 



I 



1 - e-l 3 
1 



-T+l) 



In 



(34) 
(35) 



Confirming the solution of a multiple disease Penna 
model through simulation presents considerable prob- 
lems. The solution presented here is not unique and there 
is no reason for a given simulation to settle on this par- 
ticular steady-state distribution. Within any finite pop- 
ulation it has been demonstrated |23| that if individuals 
are arbitrarily labelled, then after sufficient time, the en- 
tire population will have descended from individuals with 
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FIG. 5: Lifespan distribution for i max = 30, T = 4. 



just one of these arbitrary labels. The ansatz used treats 
the non-terminal mutations as arbitrary labels so the as- 
sumption that they are evenly distributed will be upset 
by finite sized population dynamics. 



F. The solution of a multiple disease Penna model 
with external death rate and reproductive threshold 

ages 

As in the single mutation case, it is a relatively simple 
matter to adapt the multiple disease solution to incor- 
porate an external death rate and birth thresholds. An 
upper birth cut-off is not considered as this will only serve 
to artificially lower l max . We do not spend any time on 
the derivation here but present the results for those who 
may consider them of particular interest. 



Where nil) 



\-o l 



*r(0 
*r(0 



n(0, 1) and X 7 (Z) is given by 



for I < Ib 
for In < I 



(36) 



n(l) C l T _ x I -a 1 

eW- T +V-b Xj (l) 



(37) 



x 



eW+i-T+i) _ bx ^i + i)( e /3 +T-1) 



In 



'ill ;") v < 



In a 

/3(imax-T + l) 



(38) 
(39) 



G. A Continuum Penna model 

The simple Penna model can be reformulated so that 
rather than considering discrete time steps, time is 



treated as a continuous variable t. As before any sub- 
population is characterised by its age a and string length 
I that are no longer constrained to be integers. The birth 
rate b is now the probability of an organism reproducing 
in unit time. Likewise, the mutation rate, m, is the prob- 
ability of a mutation occurring in unit string length. The 
steady-state equation for this model is 



n(0,0 = be' 131 / n(x,l)dx (40) 
Jo 

p oo poo 

+ mbe- f31 / dx I dl'n(x,l'). 



As in the discrete case the sum over ages, now an inte- 
gral, can be evaluated employing the fact that n(0, 1) is 
constant over time. The integral over ages of n(x, I) can 
be written as n(l), where n(l) = I x n(0, 1), giving 

Q = be- pl n{l) + - T ^-+mbe- 131 f n(l')dl'. (41) 
' J i 

The steady-state equation can be re-written in the form 
below and the integrals evaluated numerically. 

n{l + x) _ l + x e^-bl J{+* -f^di' 

n(l) " I eW+*) -b(l + x) ' V ' 

As in the discrete Penna model, in the steady-state, 
^Gmax) must be self sustaining and all sub-populations 
with a lower gentic lifespan partly reliant on mutation. 
As I is no longer constrained to being an integer the con- 
ditions become 



1 



(43) 

Gmax - S)be-W™ < 1 (44) 
Where 5 is arbitrarily small. These conditions give 

1 



Lax ^ 



b = 



1 



(45) 
(46) 



H. 



A continuum Penna model with multiple 

diseases 



A multiple disease Penna model can also be reformu- 
lated into a continuum case. The mutations are consid- 
ered to be 5 functions an d as in the discrete case, an 
organism dies once it has accumulated T diseases. 

Recalling the first-order steady-state equation from the 
discrete multiple mutation model 



= be-M- T+ Vn T (l) 



ut(1) 



(47) 



Tmbe-^ l - T+1 ^n T (l'), 



v>i 
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FIG. 6: A plot of genetic lifespan distribution for a discrete 
( x ) and continuum Penna model with Z max = 30. 



In the continuum model, where mutations no longer take 
up a finite length on the string, this will become 



= be- pl nr{l) 



n T (l) 



(48) 



+ Tmbe- pl [ n T (l')dl'. 
Jv>i 



An integral equation can be derived from this using 
a similar approach as in the single disease continuum 
case. To evaluate n(lr) all possible configurations of non- 
deleterious mutations must be summed over. 
In the discrete Penna model this leads to 



n(i) = C l T _ inT {l), 



(49) 



as each mutation takes up one site on the string. When 
positioning 5 functions, there is no possibility of muta- 
tions overlapping as they are of infinitesimal size on the 
string. In the continuous case, n(Z) is given by 



n(l) = V^n T {l) 
ii(It + x) is given by 

n{l + x) (l + x) T - 1 l + x 



n(l) 



IT-1 

e? 1 - bl 



(50) 



(51) 



J3(l+x) _ b Q + x j 



fl + x rnbTl' dl , 



The steady-state conditions are unchanged from the sin- 
gle disease continuous Penna model as T > 1 will only 
affect contributions from mutation that have no effect on 
the sub-population with string length Z max . The steady- 
state conditions are 



1 

< - 



b = 



1 



(52) 
(53) 



III. A PENNA MODEL WITH ARBITRARY 
BIRTH AND SURVIVABILITY FUNCTIONS 

The methods employed to solve the variety of Penna 
models discussed so far can be generalized to any Penna 
model where the birth and survival functions are func- 
tions of the organisms' age and genetic lifespan. The sur- 
vival function is defined so that the number of organisms 
of age x is given by n(x,l) — f s (x,l)n(0,l). The birth 
function b(x, I) gives the average number of offspring pro- 
duced per time step by an organism of age x and string 
length I. the steady-state equation can be written as 



n(0,0 = e~ fil ]T b(x, l)n(x, I) 



(54) 



x=0 

oo oo 



mi - b(x, l')n(x, I'). 

x=0 V>1 



Defining n(l) = J2^L n(x,l), where n(l) is the number 
of organisms with string length / regardless of their age, 
and employing and L{1) where 

oo 

L{1) = £/«(*,0. (55) 

oo 

b X (l) = *)/.(*> 0. (56) 



x=0 



the steady-state equation can be written as 



^-'W) n(,) -W) +mbe ~ S %W) nn (57) 

L(l) is the expected lifespan of an organism with string 
length I and b\(l) is the expected number of offspring 
from an organism of string length I throughout its life, b 
has been chosen so that throughout this paper the highest 
birth rate in any time step is b. Using the steady-state 
equation a general recursion relation can be generated 



X(0 nn 



n(l + 1) _ L(l + 1) [eP l - b X (l)} 



i(0 ' L(l) [e^ 1 ) - b X (l + l)e~P] 



(58) 



The conditions imposed on Z max to ensure that n(Z max ) 
is self sustaining and that the population remains finite 
are 



Thus 



M^max)e^ max = 1, 

&X(U« - l)e-^— - 1 ) < 1. 



X^rnax 1) p 



X(^max) 



(59) 
(60) 



(61) 
(62) 
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The survival function f s (x,l) should for physical rea- 
sons be a monotonically declining function. There are no 
constraints on the birth function other than that it must 
be positive. Through suitable choices of birth and sur- 
vivability functions, the Penna model can be adapted to 
model a wider variety of real-life behaviour. In previous 
work we have demonstrated that a subtle modifi- 
cation to the simple Penna model, namely replacing the 
step-function survivability with a Fermi function, is ca- 
pable of producing a mortality plateau. This particular 
modification does not change the steady state equation 
from that of the Simple Penna model, though this would 
not pose any problem to the general solution presented 
here. 

Incorporating the methods used to solve continuous 
and T > 1 models in this paper, will allow for application 
of this general method to arbitrary T > 1 and continuum 
cases. 



IV. EXTRACTING MORTALITY DATA FROM 
GENETIC LIFESPAN DISTRIBUTIONS 

The solutions presented so far have derived relation- 
ships for distribution of genetic lifespans of a population. 
From genetic lifespan distributions mortality behaviour 
can be evaluated. Recall that a distribution was deter- 
mined for n(l), that is the number of organisms with 
genetic lifespan I. The required quantitity for evaluating 
mortality rates is n(0, 1) and is given by 7^17 ■ 

The number of organisms of genetic lifespan I dying 



between age x and x + 1 is given by 

n(x, I) - n(x + 1,1) = n(0, l)[f.(x, I) - f.(x + 1, 1)]. (63) 

The number of organisms dying at age x is this quantity 
summed over all genetic lifespans. The fraction of the 
total population dying at age x is defined as the time- 
step mortality rate A4(x) 



M( Ein(0,l)[fs(xj)-fs(x+l,l)} 
1{X> £„ n((y') W) 



(64) 



V. CONCLUSION 

We have presented an analytic solution to the steady- 
state Penna model capable of dealing with arbitrary sur- 
vival and birth functions. We hope this will stimulate 
further modifications to the Penna model, where suit- 
able choices of birth and survival functions will allow an 
adapted Penna model to encompass and explain more 
observed phenomena within age structured populations. 
In our future work we will study the dynamics of the 
Penna model and consider more sophisticated complica- 
tions such as the introduction of a positive mutation rate. 

The authors have been assisted and inspired by Pro- 
fessor M. E. Cates, Professor L. Demetrius, Dr. R. Farr, 
Dr. M. J. Rutter, and Professor D. Stauffer who are grate- 
fully acknowledged. 
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